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Classical Hamiltonian spin systems are continuous dynamical systems on 
the symplectic phase space (5 2 ) n . In this paper we investigate the underly¬ 
ing geometry of a time discretization scheme for classical Hamiltonian spin 
systems called the spherical midpoint method. As it turns out, this method 
displays a range of interesting geometrical features, that yield insights and 
sets out general strategies for geometric time discretizations of Hamiltonian 
systems on non-canonical symplectic manifolds. In particular, our study 
provides two new, completely geometric proofs that the discrete-time spin 
systems obtained by the spherical midpoint method preserve symplecticity. 

The study follows two paths. First, we introduce an extended version of 
the Hopf fibration to show that the spherical midpoint method can be seen 
as originating from the classical midpoint method on T*R 2n for a collective 
Hamiltonian. Symplecticity is then a direct, geometric consequence. Second, 
we propose a new discretization scheme on Riemannian manifolds called the 
Riemannian midpoint method. We determine its properties with respect to 
isometries and Riemannian submersions and, as a special case, we show that 
the spherical midpoint method is of this type for a non-Euclidean metric. In 
combination with Kahler geometry, this provides another geometric proof of 
symplecticity. 
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1 Introduction 


A well-known integrable PDE is the continuous classical Heisenberg equation of ferro¬ 
magnetics 

w = w x w", w: S 1 —> 5 2 , 

where we represent elements in S 2 as unit vectors in R 3 . Spatial discretization of this 
equation by 

, w(s — As) — 2w(s) + w(s + As) 
w (»)=»-- 

leads to the classical Heisenberg spin chain 

Wi = WiX (Wi -1 + tOi+l), Wi E S 2 , w 0 = w n . 


More generally, a classical Hamiltonian spin system is of the form 


Wi = Wi X 


d H 

dwi 1 


Wi E S 2 , i = 1,..., n, 


( 1 ) 


for some Hamiltonian H : ( S 2 ) n —> R. In addition to ferromagnetics, examples include 
the free rigid body and the motion of n point vortices on the sphere. A key property 
of the flow of (1) is preservation of the symplectic structure of ( S 2 ) n . In the literature 
on spin systems it is common to think of the Hamiltonian as a function H : R 3n —?> R. 
Notice, however, that only its restriction to ( S 2 ) n affects the dynamics. 

We are interested in symplectic, discrete-time versions of the equations (1). Towards 
this goal, a fruitful approach is to regard the two-sphere S 2 as a coadjoint orbit of the 
Lie-Poisson manifold so(3)* corresponding to the Lie algebra so(3) of skew-symmetric 
matrices (for details on Lie-Poisson manifolds, see [7] and references therein). Then 
one can use variational discretizations, as those developed by Moser and Veselov [17] 
for some classical integrable systems, particularly the free rigid body (see also [15, 13] 
for the extension to arbitrary Lie-Poisson manifolds and Hamiltonians). This discrete 
Moser-Veselov (DMV) algorithm is formulated as an SO(3)-symmetric symplectic map 
on the phase space T*SO(3). The symmetry implies that the flow descends to a flow 
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on T*SO(3)/SO(3) ~ so(3)*, but an explicit representation on so(3)*, without using 
auxiliary variables in SO(3), is not available. Therefore, an open problem has been to 
find a minimal-coordinate symplectic discretization of the equations (1). A solution is 
provided by the spherical midpoint method , first communicated in [12]. This method is 
given by the map (w i,.. ., w n ) H > (Wi,..., W n ) implicitly defined by 

Wj - Wj _ Wi + \Vj 8H / wi+Wi w n + W n \ 

At \wi + Wi\ X dwi + Wi \’ ” ’ ’ \w n + W n \J' 

A direct proof of its symplecticity is given in [11], where also several examples for specific 
Hamiltonians are given. The proof in [11] is via a lengthy direct calculation that is not 
too enlightening. In this paper we carry out an in-depth geometric investigation of the 
method (2) and the corresponding discrete-time spin systems. 

Although there has been extensive interest in Lie-Poisson integration, and in the 
associated discrete mechanics [22, 8, 20, 14, 6, 4, 21], all previous methods have been 
closely related to the classical generating functions defined on symplectic vector spaces, 
and all use extra variables. The map defined by (2) uses no extra variables and is 
in some sense the first generalization of the Poincare generating function [19, vol. Ill, 
§319] (corresponding to the classical midpoint method) to a noncanonical, nonlinear 
phase space. Earlier discrete Lie-Poisson mechanics has also led to interesting new 
discrete integrable systems, including the hugely influential Moser-Veselov system [1, 2, 
3]. The method (2) applied to the free rigid body leads to an integrable mapping of 
an apparently new type. These considerations motivate our study of the geometry and 
discrete mechanics associated with the spherical midpoint method. 

First, in § 2, we study the symplectic geometry. In particular, we show that the 
spherical midpoint method can be interpreted as a collective symplectic integrator, such 
as developed in [9, 10]. We establish this connection by an extension of the classical Hopf 
fibration. In addition to geometric insights, the connection to collective integrators also 
establishes an independent geometric proof of symplecticity, completely different from 
the direct proof in [11], 

Second, in § 3, we study the Riemannian geometry. The classical midpoint method 
evaluates the vector held at the midpoint of a straight line joining the start and end¬ 
ing points. This suggests a generalization to Riemannian manifolds, and a Riemannian 
midpoint method, that appears to be new. We introduce this method and establish some 
of its basic properties, including equivariance with respect to the isometry group of the 
manifold and natural behaviour with respect to Riemannian submersions. Perhaps coun¬ 
terintuitively, the Riemannian midpoint method for the standard Riemannian structure 
on S 2 is not symplectic. Nevertheless, there is a Riemannian metric for which the corre¬ 
sponding Riemannian midpoint method is exactly the spherical midpoint method. We 
arrive at this result by examining how the Kahler structure of the space of quaternions 
relates to the extended Hopf map. This provides another way to view the spherical 
midpoint method, and yet another proof of symplecticity, based on Kahler geometry. 

We use the following notation. X(M) denotes the space of smooth vector fields on a 
manifold M. If M is a Poisson manifold, and H E C°°(M) is a smooth function on M, 
then the corresponding Hamiltonian vector field is denoted Xh- Let us also recollect the 
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concept of intertwining. To this extent, let M and N be two manifolds, and consider a 
differentiable map f:N—t M. We say that f intertwines X E £(M) and Y E X(N) if 
X o f = Tf o Y. (Some authors prefer to say that X and Y are / -related.) Likewise, we 
say that / intertwines a function $: M — > M and a function T: N —> N if 

o / = / o T. 

Finally, the Euclidean length of a vector w E R d is denoted \w\. If w E R 3n ~ (R 3 )", 
then Wi denotes the Lth component in R 3 . 

We continue this section with a concise presentation of the spherical midpoint method. 


1.1 Spherical midpoint method 

Here we review some background on the spherical midpoint method (2). All the results 
in this section are also available in [11]. 

A key point in this paper is the relation between the spherical midpoint method and 
the classical midpoint method on vector spaces. We recall its definition. 


Definition 1.1. Let A be a vector field defined on an open subset of a vector space. 
The classical midpoint method for X is the mapping z i—^ Z defined by 


= Y ( Z + z \ 

At V 2 /’ 


( 3 ) 


where At > 0 is the time step. 

The vector field Xh given by the right-hand side of (1) is defined on ( S 2 ) n (since 
the Hamiltonian H is defined on (S' 2 )"). In order to relate the spherical to the classical 
midpoint method we need to extend the vector field Xu to (R 3 \{0}) n . For this, we 
define a projection map p by 



The basic observation is then that the spherical midpoint method (2) can be written 


W-w 
At 



W + w 
2 


( 4 ) 


Comparing with (3), we see that the spherical midpoint method is the classical midpoint 
method applied to Xu o p. This observation is the starting point for our developments. 

What is then characteristic for vector fields of the form £o p with £ E £((S 2 )")? This 
question leads us to the next cornerstone in the paper. 

Definition 1.2. The ray through a point w E (R 3 \{0}) n is the subset 


{(Altai,..., X n w n )- A E R+}. 
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The set of all rays is in one-to-one relation with ( S 2 ) n . Note that the vector held X = 
£ o p is constant on rays. The property of being constant on rays is passed on from 
Hamiltonian functions to Hamiltonian vector fields. 


Lemma 1.3 ([11]). If a Hamiltonian function H on (R 3 \{0}) n is constant on rays, then 
so is its Hamiltonian vector field Xh, defined by 

k=X K 

The implication is that we may replace S 2 with the manifold of rays, and Hamiltonian 
functions on ( S 2 ) n with Hamiltonian functions on R 3 \{0} that are constant on rays. 
In this representation, the spherical midpoint method becomes the classical midpoint 
method, as we have seen. 

R 3n and (R 3 \{0}) n are Poisson manifolds with the Poisson bracket 


{F,G}(w) = J2i 

k= 1 


dF(w) 

dw k 


x 


d G{w) 
dw k 


) • w k . 


( 6 ) 


This is the canonical Lie-Poisson structure of (so(3)*) Tl , or (su(2)*) n , obtained by iden¬ 
tifying so(3)* ~ R 3 , or su(2)* ~ R 3 . For details, see [7, § 10.7] or [9]. 

The flow of a Hamiltonian vector held Xh on R 3n , denoted exp (A#), preserves the 
Lie-Poisson structure, i.e., 

{F o exp (X h ),G o exp (X H )} = {F, G} o exp(X H ). VF, G G C°°(R 3n ). 


The how exp(Xn) also preserves the coadjoint orbits [7, § 14], given by 

S 2 Xl x • • ■ x S 2 Xn c R 3n , A = (Ai,..., An) G (R + ) n , 

where S\ denotes the 2-sphere in R 3 of radius A. A Lie-Poisson integrator for Xh is a 
method that, like the exact how, preserves the Lie-Poisson structure and the coadjoint 
orbits. 

It is possible to extend the spherical midpoint method so that it encompasses all 
non-singular coadjoint orbits, instead of only the one with radius one. Define the map 
T(w,W) by 


r: 


V \w\ + W\\ 


VT 


\W n \(w n + W n ) 

W n T n | 


We then have the following dehnition. 

Definition 1.4. Let X be a vector held on an open subset of R 3n . The extended 
spherical midpoint method for X is the discrete-time system w t—> W dehned by 


W-w 
A t 


X(T(w,W)). 


( 7 ) 
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A consequence of the geometric investigations in § 2 is that the method (7) is a Lie— 
Poisson integrator, directly related to the classical midpoint method on T* R 2n through 
the concept of collective symplectic integrators. A consequence of the geometric investi¬ 
gations in § 3 is that the method (7) is a Riemannian midpoint method with respect to 
a non-standard metric on (R 3 \{0}) n . 

2 Symplectic and Poisson geometry 

In this section we show that the extended spherical midpoint method and the classical 
midpoint method are coupled through collective symplectic integrators [9, 10]. 

We use quaternions as it makes the calculations more transparent; the field of quater¬ 
nions is denoted H. We apply the convention that product sets C n and H n inherit the 
componentwise operations of the underlying field. For instance, if z = (z\, Z 2 ) 6 C 2 , 
then 

z 3 = (zf,zi). (8) 

All operations are defined in the same manner. 

2.1 Intertwining by the double covering map 

We first consider intertwining in the double covering case. We define the double covering 
map 


vo: C n —»• C n , z^z 2 (9) 

following the convention in (8). We use the notation C* := C\{0} and C” := (C*)" - . 

Lemma 2.1. Let X, Y e 3£(C*) and let <h(AfX) and denote the classical 

midpoint method (3) on C n for X and Y respectively. Assume that: 

1. X(Xz) = X(z) for all A <E R£, i.e., X is constant on rays. 

2. Y is tangent to the tori in Cf, i.e., Y{z)/z is imaginary for all z e C”. 

3. vo intertwines X and Y 

Then vo intertwines <h(AfA) and $(AfY). 

Proof. The proof is illustrated in Figure 1. Consider two points z and Z, solutions of 
one step of the classical midpoint method for Y. The midpoint is z := (z + Z)/ 2, so 

Z-z = A tY(z). 

The assumption that vo intertwines X and Y is 

X(z 2 ) = 2 zY(z). 
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Figure 1: An illustration of Lemma 2.1. The midpoint z maps to the same ray as the 
midpoint z 2 := (z 2 + Z 2 )/2. 

We have 

Z 2 - z 2 = 2z(Z - z) 

= 25y(5) (10) 

= A tX{z 2 ). 

Consider the general identity 

z 2 + Z 2 = i((Z-z) 2 + (z + Z) 2 ). 

Without loss of generality, we assume that z G R n . We assumed that Y was tangent 
to circles, so Y(z) has only imaginary components, so the same holds for Z — z, which 
implies that (Z — z) 2 only has real components. We therefore obtain that z 2 + Z 2 is in 
R n . Since X is constant on the rays, and since R" is a ray in C”, we get 

X((z 2 + Z 2 )/2) = X(z 2 ). 

From (10) we now have 

Z 2 - z 2 = AtX ((z 2 + Z 2 )/2). 

This proves the result. □ 




2.2 Intertwining by the extended Hopf map 

Consider the map 

vr: H” —» H n 

z i —y -zkz 
4 
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Again, we follow the convention of (8) and all the operations are defined componentwise. 
Note that the image of ir has no real part. Let us define the three dimensional subspace 
of pure imaginary quaternions by 

V := span{ i, j, k } 

If we identify V with R 3 , we can regard tt as a map 

7r: H n —)• R 3n . (11) 

When n = 1 this is the extended Hopf map, essential in the construction of collective 
Lie-Poisson integrators on R 3 [9]. 

V n is naturally endowed with the Lie-Poisson structure (6), inherited from R 3n . We 
use the notation V* ■= W\{0} and VJ 1 := (14) n . 

Lemma 2.2. Let X E X(V™) and Y E X(H"), and let $>(AtX) and ^f(AtY) denote the 
classical midpoint methods on V n and H n respectively. Assume that: 

1. X(Xw) = X(w ) for all A E R” and w E VJ 1 '. 

2. Y is tangent to 3-spheres, i.e., z~ x Y(z) E V n for all z E H". 

3. Y is orthogonal to the fibres of n, i.e. kz~ 1 Y(z) E V n for all z E H™. 

4- tt intertwines X and Y. 

Then tt intertwines 3>(A tX) and 'fl(AtY). 

Proof. Consider z and Z, solution of the classical midpoint method in H n for Y. The 
midpoint is denoted by 

„ z + Z 

z := -. 

2 

Since the classical midpoint method is equivariant with respect to affine transformations, 
we may, without loss of generality, assume that z is real (i.e., all its components are 
real). At the point z, the real direction is orthogonal to the 3-spheres, and k is the fibre 
direction. Without loss of generality, we may further assume that Y(z) is proportional 
to i, i.e., Y(z) = i a and a E R". As a result, the components of z and Z belong to 
span{l,i}, which we identify with the complex plane C, so we write z, Z E C n . 

Notice that since ki = —ik, we have for z = a + ib: 

zk = z k 

When restricted to C n , the Hopf map at z E C n is thus 

7 r(z) = zkz = z 2 k 

This means that 7r(C n ) C (spanjk,j}) n . We identify (span{k,j}) n with C n by k 1 
and j 0 —i. With these identifications, the restriction of tt to C n is the double covering 
map w defined in (9), so the result follows from Lemma 2.1. □ 
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H* carries the structure of a symplectic manifold; the Hamiltonian vector field corre¬ 
sponding to F E C^H”) is 

X F (z) = VF{z)k. (12) 

The symplectic structure coincides with the canonical symplectic structure of T*R 2n 
under the identification 


rl* 3 a + 10 + JC 




7 / 7 


P 


Likewise, Vfi carries the structure of a Poisson manifold; the Hamiltonian vector field 
corresponding to H £ C°°(Vfi) is 


X H (w) = n (wVH(w)), 


(13) 


where n: H n —> V n is the projection a + ib + jc + kd i->- ib + jc + k d. Under the 
identification of V with R 3 , this Poisson structure coincides with the standard Lie- 
Poisson structure of R 3 , so equation (13) is just another way of writing equation (5). 

We now investigate what the conditions on X and Y in Lemma 2.2 mean for Hamil¬ 
tonian vector fields X = Xu and Y = Xp. It follows directly from Lemma 1.3 that 
Xu fulfils condition 1 in Lemma 2.2 if and only if H fulfils the same condition, i.e., 
H(Xw) = H(w) for all A 6 R" and w 6 The next result shows that Xp fulfils 
condition 3 in Lemma 2.2 if and only if F fulfils the same condition. 

Lemma 2.3. Let F E C°°(H0). Then Xp is orthogonal to the fibres of n, i.e., 

k z~ 1 X F (z) E V n 


if and only if 


F(Xz) = F(z), VAeR" and z E H£. 


Proof. From (12) it follows that Xp(z) is orthogonal to the fibres if and only if 


kz~ 1 XF(z)k E V n . 


This is equivalent to 

z~ l XF(z) E V n , 

which means that z~ 1 'VF(z) is pure imaginary, so XF(z) is tangential to the spheres. 
Since this is true for any z E H”, it means that F(Xz) = F(z) for A E R". □ 

Given H E C' 00 ^"') we can now construct discretizations of systems on Vfi in two 
ways: 

1. The classical midpoint method on H” for the vector field Xu on descends to a Lie- 
Poisson integrator on Vfi. The resulting discrete-time system, examined in [10], is 
an example of a collective symplectic integrator. 
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2. The extended spherical midpoint method (7) for the vector field Xu gives a 
discrete-time system on Vff. 

In general the two methods are different. They do, however, coincide for ray-constant 
Hamiltonian vector fields, which is the main result of this section. 

Theorem 2.4. Let H £ C°°(V™) be constant on rays, let denote the classical mid¬ 
point method on V™, and let T denote the classical midpoint methods on H*. Then the 
extended Hopf map it intertwines &(Xu) and 4 , ( Xuon )• That is, the map on Vf 1 induced 
by ^(X Hon ) coincides with &(Xu)- 

Proof. Appealing to Lemma 1.3, the vector field Xu is constant on rays. If H(Xw) = 
H(w ) for all A £ R™, then H on fulfils the same property, since n is homogeneous, that 
is: 


7t(A;z) = A 2 7 t(z), A £ R”. 

We can therefore use Lemma 2.3 to obtain that Xu on is orthogonal to the fibres. 

The Hopf map is a Poisson map [9], so 

T z ir ■ Xuon(z) = X h (tt(z)), ze H, 

which means that it intertwines Xu and Xu on- The result now follows from Lemma 2.2, 
since 4>(A^) coincides with the classical midpoint method applied to Xu- □ 

Corollary 2.5. The spherical midpoint method, defined by (2), is symplectic map and 
the extended spherical midpoint method, defined by (7), is a Lie-Poisson integrator. 

Proof. Follows from Theorem 2.4 and the result in [10] that collective symplectic inte¬ 
grators are symplectic. □ 

3 Riemannian and Kahler geometry 

In this section we describe the geometry of the extended spherical midpoint method 
from the viewpoint of Riemannian and Kahler geometry. More precisely, we construct a 
method on Vf 1 , stemming from a non-Euclidean metric, that coincides with the spherical 
midpoint method for Hamiltonian functions that are constant on rays. The relation 
between these two methods is established through the classical midpoint method on H”. 
Before working this out in §3.2, we develop, in §3.1, a general theory of midpoint methods 
on Riemannian and Kahler manifolds. This theory further reveals the geometry of the 
spherical midpoint method, and provides a starting point for generalizations. 

3.1 Riemannian midpoint methods 

Given a Riemannian manifold (M, g), let [0,1] B t t-> 7 g (t]W,W) £ M denote the 
geodesic curve between w and W. 
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Definition 3.1. Given a vector field X £ X(M), the Riemannian midpoint method on 
M is the discrete-time system <3? g (A tX) : w H > W defined by 


d_ 
d t 




AtX( lg (l/2,w,W)). 


(14) 


If M = and g is the Euclidean metric, then (14) coincides with the definition of 
the classical midpoint method (3). 

Riemannian midpoint methods transform naturally under change of coordinates: 


Proposition 3.2. Let M and N be two diffeomorphic manifolds, let ip: N M be a 
diffeomorphism, and let g be a Riemannian metric on M. Then 

'0“ 1 o 4> g (AtA) o if) = <h^,* g (A tip*X). 


Proof. The result follows from the definition (14) of <h g and standard change of coordi¬ 
nate formulas in differential geometry. □ 


A consequence of Proposition 3.2 is that Riemannian midpoint methods are equivari- 
ant with respect to isometric group actions: 

Proposition 3.3. Let (M, g) be a Riemannian manifold, and let G be a Lie group 
acting isometrically on M. Then the Riemannian midpoint method 4> g is equivariant 
with respect to G, i.e., 

g - 1 ° $g( AtX ) ° V’p = $g {AtlfgX) 
where r ijj g : M —> M denotes the action map of g £ G. 

Proof. The result follows from Proposition 3.2 and ip* g = g (the action is isometric). □ 

We will now discuss a generalised version of Proposition 3.2, where M and N are no 
longer diffeomorphic. 

Let 7 r: N M be a submersion from N to another manifold M (tv is smooth and its 
Jacobian matrix T z n: T Z N — » T n ^M is surjective at every z £ N.) Then ir induces 
a vertical distribution Vert by Vert z = {n £ T z N\T z ir ■ v = 0}. By construction, the 
vertical distribution is integrable, and the fibre through z G N is given by 7 r _ 1 ({ 7 r(z)}). 
If ( N , h) is Riemannian, then the orthogonal complement with respect to h is called 
the horizontal distribution and denoted Hor. Typically, the horizontal distribution is 
not integrable. The Riemannian metric h is called descending (with respect to the 
submersion 7 r) if there exists a Riemannian metric g on M such that for all z £ N 

h z (u,v) = g^( z )(T z vr • u,T z n ■ v), Vu,v £ Hor z . 

The map 7 r between the Riemannian manifolds ( N , h) and (M, g) is then called a Rie¬ 
mannian submersion. For details on Riemannian submersions, see [18, §3.5]. 

A vector field Y £ X(N) is called horizontal if Y(z) £ Hor z for all z £ N. Y is called 
descending if there exists a vector field X £ X(M) such that 7 r intertwines X and Y, 
i.e., T z ir ■ Y(z) = X(i r(z)). The following result, schematically illustrated in Figure 2, 
is a generalized version of Proposition 3.3. 
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Figure 2: An illustration of Theorem 3.4. If the vector fields Y and X are intertwined 
by the Riemannian submersion n and Y is horizontal, then the Riemannian 
midpoint methods for X and Y are also intertwined by n. 


Theorem 3.4. Let (M, g) and (IV, h) be Riemannian manifolds and n: N —> M a 
Riemannian submersion. Let Y e X(N) be horizontal, let X 6 X(M), and assume 
that tt intertwines X and Y. Then n intertwines the Riemannian midpoint method 
<h g (hX) and the Riemannian midpoint method 4>h(AtT), i.e., 

7r($ h (A tY)(z)) = $ g (AtX)(7r(z)). 


Proof. Let z and Z fulfil (14). Let w = 7 r(z) and W = 7t(Z). We need to show that 
W = <hg(A t.X)(w). The geodesic 7 h (t\z,Z) is horizontal at t = 1/2. It is therefore 
horizontal at all times [18, § 3.5]. Since horizontal geodesics on N maps to geodesics 
on M, we have that 7r(7h(i; z, Z)) = 7 g (t; w, W). By applying Tir to (14) we obtain 


(1/2, z ,Z)7T 


d. 
d t 


t= 1/2 


7h 


Z) 


T lh (i/2,Z,Z)K ■ AtT(7h(l/2;«,Z)) 


d_ 

dl 


2 7r ( 7h ( t ’ Z ’ Z ^ = AtT 7h P/Lz,Z) 7T ’ ^(7h(l/2; 2, Z)) 


d_ 

dl 


2 7r ( 7h ( t;z ’ Z )) = ALV(7r(7h(l/2 \z,Z))) 


d_ 
d t 




XtX{ lg (l/2-w,W)). 


Thus, W fulfils the equation defining 4> g (AfA'), which proves the result. 


□ 
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Next, assume (N. h,w) is a Kahler manifold. Recall the properties of a Kahler mani¬ 
fold: there is a map J: TN —>• TN called the complex structure that fulfills 

h(u, v) = l v(u, Jv) 
u(u,v) = h(Ju,v) 
h (Ju, Jv) = h (u, v) 

X F = J~ l VF , F £ C°°(N) 

where V is the gradient with respect to h. We are interested in the case when vr is both 
a Riemannian submersion and a Poisson map. 

Lemma 3.5. Let (N, h,w) be a Kahler manifold, let (M, g, {•, ■}) be a Riemannian and 
Poisson manifold, and let H £ C°°(M) be a Hamiltonian. Assume there is a Riemannian 
submersion ir: N —> M that is also a Poisson map. Then Xho-k Is horizontal if and only 
ifVH (gradient of H with respect to the metric on M) is tangent to the symplectic leaves 
of M. 

Proof. By definition, the vector field Xho-k is horizontal if and only if 

h(Xeo 7 r, v) = 0, Va eVert. (16) 

By (15) we also have 

h(X Hon ,v) = h( J~ X V (H o 7r),v) 

= h(V(LTo7r), Jv) (17) 

= u(v,X(H o 7r)). 

Combining (16) and (17), Xn on is horizontal if and only if 

ui(v, X7(H o 7r)) = 0, VaGVert. 

Expressed in words, X Fon is horizontal if and only if V {H o n) belongs to the symplectic 
complement of Vert, denoted Vert^“. Let Pre denote the distribution on N defined by 
the preimage of the tangent spaces of the symplectic leaves of M. From [5, Proposi¬ 
tion III.14.21] it follows that Pre = Vert^. It remains to show that Vi7 is tangent to 
the symplectic leaves if and only if V (H o n) G Pre. Since the metric h is descending, 
the gradients on M and N are related by 

Tir o X(H 0 7r) = VH. 

By the definition of Pre, this formula proves the result. □ 

By combining Theorem 3.4 with Lemma 3.5 we obtain the following result. 

Theorem 3.6. Let M, N, and 7r be as in Lemma 3.5. Let H £ C°°(M) fulfil the 
condition in Lemma 3.5, i.e., X7H is tangent to the symplectic leaves. Let 4> g and 
denote the Riemannian midpoint methods on M and N respectively. Then: 
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1. If &h{XHon) is symplectic, then & g (Xff) is a Poisson map. 

2. If ^h(XHon) preserves the pre-image of the symplectic leaves, then pre¬ 

serves the symplectic leaves. 

3. If G is a Lie group that acts on M and N, and 4>h and tt are equivariant with 
respect to G, then 4> g is equivariant with respect to G. 

Proof. Let and Tg be the Riemannian midpoint methods on M and N respectively. 
By Lemma 3.5, Xho-k is horizontal, since XH is tangential to the symplectic leaves. The 
vector field Xho-it is thus descending (it descends to Xh since tt is a Poisson submersion) 
and horizontal. By Theorem 3.4, n then intertwines ^g(Xn) and <l>h(-Vf/o 7 r)- 

Proof of (1): Since tt is a Poisson map and Q^Xho-k) is a symplectic map, 4> g (Xtf) is 
a Poisson map. 

Proof of (2): Since 3>h(-X#o7r) preserves the integral submanifolds of the symplectic 
complement of Vert, and since these submanifolds project to the symplectic leaves, it 
follows from the tt intertwining property that $g(V/f) preserves the symplectic leaves. 

Proof of (3): Let X £ X(M). Let Y £ X(N) be horizontal and descending to X. Let 
g £ G Then g ■ X o n = g ■ Tiro X = Tng ■ X , since tt is equivariant. Thus, g ■ Y descends 
to g ■ X. Next, using Theorem 3.4 

d’gCg • X) o TT = TT o $ h (g ■ Y) = TT o g~ l • <L h (K) • 9 = S' -1 • $ g PQ ■ g on. 

This proves the results since tt is a submersion. □ 

3.2 Riemannian structure of the spherical midpoint method 

Our objective is to show that the spherical midpoint method (7) on Vf 1 , for Hamiltonian 
vector fields Xh £ X{Vf L ) with H of the form in Lemma 3.5, is a Riemannian midpoint 
method with respect to a non-Euclidean Riemannian metric, related to the classical 
midpoint method on H™ by a Riemannian submersion in the sense of Theorem 3.4. 

Recall that the extended Hopf map (11) is a submersion n: H" —> Vff that is a Poisson 
map with respect to the Kahler structure on H" and the Poisson structure on Vff (as 
described in § 2). 

Lemma 3.7. The Kahler metric on H" is descending with respect to the extended Hopf 
map tt : H” — > Vf l . The corresponding Riemannian metric on Vff is 

n 

g w{u,v):=Y j V ^y. (18) 

1=1 1 1 

Proof. Each fibre 7r -1 ({w}) C H" is the orbit of an action of the group U(l) n on H”. 
This action is isometric with respect to the Kahler metric. That is, if h denotes the 
Kahler metric and Lq denotes the action map, then Lgb = h. It follows from [16, 
Proposition 4.3] that h is descending. Direct calculations, straightforward but lengthy, 
confirm that it descends to the metric (18). □ 
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As a specialization of Theorem 3.4 to the case M = Vfi and N = H*, we obtain 
a relation between the Riemannian midpoint method on Vff and H™. Notice that the 
Riemannian midpoint method on H" is the classical midpoint method, since the metric 
of H” is Euclidean. We denote the classical midpoint method on H" by 4> cm . 

Theorem 3.8. Let Y E X(H”) be a horizontal vector field, let X E XiVf 1 ), and assume 
that the extended Hopf map it intertwines X and Y. Further, let g denote the Rieman¬ 
nian metric (18) on Vfi. Then n intertwines the Riemannian midpoint method 4> g (/iA') 
and the classical midpoint method Q cm (hY). 

As a specialization of Lemma 3.5 to the case M = Vff and N = H™, we obtain a 
geometric formulation of Lemma 2.3. 

Lemma 3.9. Let H E C°°(Vfi). Then Xho-k Is horizontal if and only if H is constant 
on the rays. In particular, Xho P o-k Is horizontal for any H E C°°((S 2 ) n ). 

Proof. The symplectic leaves of Vfi are the coadjoint orbits of (so(3)*) n . These consists 
of 

Sr x x • • • x S 2 n = {(tui,..., w n ) E V?; \ w k\ = r k}, (19) 

for arbitrary r & E R + . Let V denote the gradient on Vfi with respect to g. It follows 
from Lemma 3.5 that Xho-k is horizontal if and only if XH is tangent to (19). From (18) 
and the direct product structure of Vfi, we see that g (u,v) = 0 for any u tangent to 
(19) if and only if v is tangent to the rays (Definition 1.2). Let R. p ( w ) C Vfi denote the 
ray through w. The condition for X^jon to be horizontal is therefore 


g w (VH(w),v) = 0, \/v E T w R p{w) 

(d H(w),v) =0, Vu E T w R p ( w ), 

which implies that H must be constant on the rays. □ 

As a specialization of Theorem 3.6 to the case M = Vfi and N = H”, we recover 
again that the Riemannian midpoint method T g on Vfi is a Poisson integrator. 

Proposition 3.10. The Riemannian midpoint method <L g on Vff, applied to Hamiltonian 
vector fields, is a Poisson integrator. 

The final result in this section connects the extended spherical midpoint method (7) 
and the Riemannian midpoint method <h g . The two methods are different, but they 
coincide for ray-constant Hamiltonian vector fields, as a consequence of Theorem 2.4 
and Theorem 3.8. 


Proposition 3.11. Let T denote the extended spherical midpoint method (7) on Vfi 
and Tg the Riemannian midpoint method with respect to the metric g in (18). Let 
H E C^iVf 1 ) be constant on rays. Then T(A tXn) = 4> g (A tX H ). 
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Proposition 3.10 and Proposition 3.11 provide another independent proof that the 
spherical midpoint method is a symplectic discretization, this time based on Kahler 
geometry and Lemma 3.5. This proof is interesting because it allows for generalizations 
to other Kahler manifolds. 

It is a remarkable consequence of Proposition 3.11 that the non-Euclidean induced 
metric g, given by (18), has become redundant in the case when the Hamiltonian H is 
constant on rays; from (4) we see that the Riemannian midpoint method & g (AtX}{) can 
be expressed solely in terms of the classical midpoint method on V™ for such Hamilto¬ 
nians. 
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